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ABSTRACT 

It is pointed out that the meson mixing matrix elements usually considered respon- 
sible for the bulk of the observed few-body charge symmetry breaking are naturally 
g 2 -dependent in QCD. For n° - rj mixing, using the usual representation of the pseu- 
doscalar fields, the leading q 2 dependence can be explicitly calculated using chiral 
perturbation theory to one loop, the result being a significant decrease in the magni- 
tude of the matrix element in going from timelike to spacelike values of q 2 . Since it 
is the latter range of q 2 which is relevant to NN scattering and the few-body bound 
state, this result calls into serious question the standard treatment of few-body charge 
symmetry breaking contributions associated with n° - rj and p — u mixing. 
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In the standard model, charge symmetry (CS) is broken not only by electromag- 
netism (EM) but also, through the inequality of up and down quark masses, by the 
strong interactions. This strong breaking leads to important non-Coulombic contri- 
butions to CS violating observables such as the difference of low energy nn and pp 
scattering amplitudes and the A = 3 binding energy difference, ABE. 1 - 3 It is gen- 
erally believed that these non-Coulombic contributions can be accounted for in the 
meson exchange picture by isoscalar-isovector mixing (predominantly p-u>, but also 
7T° — 77) in the intermediate meson propagators of one-boson exchange graphs. 4 " 6 In the 
analysis of Refs 4-6, the p - u matrix element is obtained from experimental data on 
e + er — > 7r + 7r" in the region of the p — w interference shoulder, while the TT°-r) matrix el- 
ement is obtained by a pole model analysis of 77, rj' decays. The latter analysis (which 
has some additional shortcomings, to be discussed below) is feasible only because the 
mixing matrix elements are assumed to be ^-independent. This assumption is also 
implicit in the treatment of the A = 3 bound state and NN scattering, where time- 
like matrix elements, extracted as above, are used unchanged in the spacelike region. 
From the perspective of QCD, however, this assumption is clearly incorrect. Indeed, 
the basic hypothesis underlying the meson exchange picture is that, at low energies, 
QCD reduces to an effective hadronic theory involving mesons, nucleons, deltas (etc.). 
We know, however, that, as an effective low energy theory, such a hadronic theory 
must incorporate all possible terms, involving the low energy composite fields, which 
are not explicity forbidden by the symmetries of the underlying theory (in this case, 
QCD). In particular, there will be terms higher order in derivatives, involving also 
the quark mass matrix, which will produce a ^-dependence to all isoscalar-isovector 
mixing matrix elements. How important such effects will be for a given CS-breaking 
(CSB) observable cannot be determined without either generating the effective the- 
ory from QCD or constraining the strengths of the relevant terms in the effective La- 
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grangian from experiment. Some information, however, already suggests that these 
effects are unlikely to be negligible. First, we know that, in the case of effective chiral 
Lagrangians, higher derivative terms typically occur scaled by a dimensionful factor 
A x ~ Anf n ~ lGeV. 7 ~ w If this is true in general, one would expect significant changes 
in going from, eg., experimental data for p-uj at q 2 ~ m 2 to q 2 < 0. (We will consider 
the case of n° - -q mixing, using the chiral Lagrangian, explicitly below and see that 
this expectation is borne out.) Second, a recent model calculation 11 , in which the 
^-dependence of p - u> mixing arising from an intermediate (free) quark loop through 
the u — d constituent quark mass difference is evaluated using a monopole ansatz for 
the meson-quark-antiquark vertex function, finds a very strong q 2 dependence of the 
p - to element of the inverse vector meson propagator matrix. Such a q 2 dependence 
would translate into a significant q 2 dependence of the phenomenological p - u mixing 
matrix element (this aspect of the calculation is not treated in Ref 11-see Ref 12 for a 
detailed discussion). This result is qualitatively compatible with the discussion above 
and, despite possible reservations about its potential model-dependence, makes clear 
the dangers inherent in the neglect of ^-dependence of the usual treatment. 

One would, of course, ideally, like to reduce the model dependence of the p — oj 
result, since it is p - w mixing which produces the dominant non-Coulombic CSB 
contribution to few-body observables (in the standard treatment). Unfortunately, 
one is hampered, in this regard, by a lack of information about higher order terms in 
the (putative) effective low-energy Lagrangian involving the vector meson fields. If we 
turn our attention to the next-to-leading contribution, associated with tt° - r\ mixing, 
however, we are in much better shape, since the effective low-energy Lagrangian 
relevant to the pseudoscalar sector is well-known, to order q 4 , as a result of the 
work of Gasser and Leutwyler 8-10 . By computing n° - r\ mixing to one-loop in this 
framework, one obtains the leading ^-dependence in a straightforward, and reliable, 
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fashion. Since the effective Lagrangian is not constrained beyond order q 4 , one cannot 
determine the higher order ^-dependent contributions, which restriction limits the 
validity of the results to the range \q 2 \ < 0(m 2 ), by the usual power counting arguments 
of the low-energy expansion. 

Let us consider, therefore, n° - -q mixing. Recall that the relevant terms in the 
effective chiral Lagrangian, to order q 4 , are given by 8 

L eff = - A f 2 Tr(d^d^) + \fTr\pLM{H + £+)] + U [Tr(d^d^)] 2 

+ L 2 Tr{d^d v ^)Tr{d^d v ^) + L 3 Tr(d^d^Y,d v E^d v E) 

+ L 4 Tr{d ll Xdi J -tf)Tr[2iJ,M(Y: + E+)] + L b Tr [2/i(MS + ^M)d^d^] 

+ L 6 [Tr[2 M M(S + £t)]] 2 + Lj [T r [2 M M(£ - St)]] 2 

+ J L 8 T 1 r[4 M 2 (MSMIl + MI] t Mi; t )] (1) 
where \i is a mass scale related to the value of the quark condensate 8 , 

S = exp(i\.Tr/f), (2) 

with A the usual SC/^) Gell-Mann matrices, 7? the pseudoscalar octet fields and / a 
dimensionful constant, equal to f n in leading order. In Eq (1), M is the current quark 
mass matrix and we have set the external fields to zero. The (scale-dependent) co- 
effieients Li, ■ ■ ■ , L$ are as defined in Ref 8. In what follows we ignore EM contributions 
to tt° --q mixing, since these are known to vanish in the chiral limit 13 and estimates for 
the effect of the departure from the chiral limit, based on the size of the relevant EM 
chiral logarithms 14 , show them to be negligible relative to the strong contributions 
resulting from L eff in Eq (1). From Eq (1) one may compute the inverse n 3 -7r 8 prop- 
agator matrix to one-loop as a function of q 2 (where q is the external four-momentum 
and 7T 3 , 7r 8 are the unmixed octet fields that would be identical to the tt°, rj fields, 
respectively, in the limit of exact isospin symmetry). Diagonalizing this matrix one 
obtains the <j 2 -dependent mixing angle for the physical propagating modes 
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(3) 



where m u ,m d ,m s are the u, rf, s current quark masses, m = (m u +md)/2, L[ are the (scale- 
dependent) renormalized low-energy constants of Ref 8 and l P = m 2 P £n{m 2 P l 'fi%), with P 
any of the pseudoscalar mesons and \i Q the renormalization scale. The q 2 -dependence 
of 9 is unavoidable at the one-loop level. Following Ref 8 one may rewrite Eq (3) 
using one-loop results for other physical observables, obtaining 



8(q 



■> V3(m 2 K O - m 2 K+ ) QC D 
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(m 2 K - ml) 
1 / - 2 



1 + Agmo 



+ (w)( 1 + (S^b) f "<'"-^»)] <4) 

where A GM o is the (electromagnetically corrected) Gell-Mann-Okubo discrepancy 

A GMO = (4m 2 K - ml - 3m 2 ) / \m% - m\). (5) 

In Eq (4), (Am 2 K ) QC D is the contribution to the K° - K + mass-squared splitting re- 
sulting from m u ^ m d . This can be extracted from the physical splitting if one knows 
the corresponding EM contribution. The latter may be considered to be known if we 
accept Dashen's theorem 13 for the EM self-energies, 

(Am 2 K ) E M = {Am 2 ) E M (6) 

since the quark-mass-difference contribution to the ir splitting is 0[(m d - m u ) 2 } and 
known to be small 15 . Dashen's theorem is, however, valid only to lowest order in the 
EM chiral expansion. Constraints on the size of potential violations are discussed in 
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Ref 16. The resulting uncertainties in the size of (Am 2 K ) E M lead to a corresponding 
uncertainty of order 20% in the overall scale of the result (4). Apart from this overall 
scale uncertainty, the result (4) should provide an accurate representation of 9(q 2 ), 
subject to the restrictions of the validity of the low energy expansion, i.e. for values 

Of |<z 2 | < 0(m 2 ). 

Before proceeding, we must point out one important feature of the calculation 
which bears crucially on the proper interpretation of its results. Recall that L e // 
is only one of a family of possible effective Lagrangians, all others being related 
to that in Eqn (1) by redefinitions of the octet fields which leave the one particle 
singularities, and hence the S-matrix elements, of the theory unchanged 17 ) . The 
propagator matrix, and hence the mixing angle, 9(q 2 ), of course, depends on the 
particular definition of the octet fields employed: as such it is not a physical quantity. 
Only when combined with a full low-energy Lagrangian including the nucleon fields 
could one obtain results for physical processes (such as, eg., NN scattering) which 
were independent of the particular choice of the meson fields. In the context of few- 
body meson-exchange model calculations we are not, at present, able to carry out this 
program and, in fact, because of the insertion of phenomeno logical form factors for the 
meson-nucleon vertices, it is actually ambiguous as to exactly what definition of the 
meson fields is being used (in the sense above). The relevant question, in this setting, 
is then whether or not the q 2 dependence of the mixing (which must be present, on 
general grounds, for any choice of fields) is typically negligible or not. If we use the 
definition of the fields which leads to L eff of the form given in Eqn (1), we find that 
q 2 dependence of the mixing is not negligible in few-body systems, and hence that 
the assumption of q 2 independence of the standard treatments is not a reasonable 
one to make. Unless, however, one formulates the entire scattering (or bound state) 
problem under consideration in a well-defined effective field theoretic manner, where 
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one can determine that the same choice for the meson fields is being used for the 
meson-nucleon vertices as in the purely mesonic sector, one is not guaranteed that 
simply plugging the result of Eqn (7) into an existing few-body code will lead to a 
reliable estimate of the CSB associated with ir° - rj mixing. This statement is, of 
course, equally true of the usual treatment, with, however, the added caveat that the 
usual treatment suffers from the absence of a demonstration that there even exists 
a choice of the meson fields for which the mixing is ^-independent. The result of 
Eqn (7) is, of course, perfectly reliable for the choice of meson fields implied by 
Eqn (1); the problem is that possible additional q 2 dependent CSB contributions 
associated with the meson-nucleon couplings cannot be consistently evaluated owing 
to the phenomenological treatment of the meson-nucleon vertices. For this reason, 
quantities such as the n° - r\ and p - to matrix elements used in conventional few-body 
CSB calculations are actually ambiguous beyond leading order. This means that 
the appropriate (and most conservative) interpretation of the present calculation is 
that the difference between the results of using the g 2 -dependent mixing of Eqn (7) 
and the usual ^-independent mixing represents a lower bound for the uncontrollable 
theoretical uncertainties associated with the phenomenological modelling of the input 
meson-nucleon vertices. The comments which follow are to viewed in this light. 

From the results above we obtain, for the choice of fields implied by Eqn (1), two 
pieces of information: first, the size of the 7r°-?7 mixing and, second, its q 2 dependence. 

Let us first consider the magnitude. The result (4) leads to a phenomenological 
mixing matrix element given by 

m 2 s (q 2 ) = -(m 2 l - m 2 7r )9(q 2 ). (7) 

As mentioned above, there is an overall uncertainty in the result (7) as a consequence 
of the unknown size of possible violations of Dashen's theorem. Since this uncertainty 
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also enters the determination of A GM o the effect is not quite linear in (Aiti^qcd, but 
it is very nearly so. EM and strong contributions to Am 2 K have the opposite sign so 
that EM contributions larger than those corresponding to Dashen's theorem lead to 
larger values of (Am 2 K ) QC D- Let us write 

mj 8 = -[a a + (q 2 /m 2 v )a 1 ]. (8) 

The quantities a Q ,ai are listed in Table 1 for a range of values of (Am K ) E M between 
1.3MeV (arising from Dashen's theorem) and 2.8MeV (which results from attempting 
to saturate the Cottingham formula for the kaon system with K and K* intermediate 
states 18 ). 

The values of Table 1 should be considered to supercede those of earlier work 
on the subject for the following reasons. The tt° - rj matrix element of Ref 6c was 
extracted by an analysis of r/, rf decays in which these decays are taken to be mediated 
by rj, rj' and tt poles, intermediate 7r°?7, tt t]' vertices being taken to be ^-independent. 
There are two problems with this analysis. First, as we have seen, the ir° - rj vertex 
(and hence, presumably, also the ir° - rj' vertex) is not ^-independent. Second, as 
one discovers from Ref 10, which treats rj — > 3n to one loop in chiral perturbation 
theory, a large portion of the 55% increase of the r\ — > in amplitude in going from 
tree level to one-loop level is attributable to enhancements associated with final state 
interactions produced by the presence of s-wave / = -kit pairs. As an s-wave tttt effect 
this enhancement can clearly not be represented by an intermediate pseudoscalar pole 
(the effect of the rj' in L eff is contained entirely in the low-energy constant L 7 8 ~>); it 
has, therefore, not been included in the pole model analysis. Its inclusion would 
presumably lower the extracted values of the n° - rj and tt° - r{ matrix elements, even 
if one could treat them as g 2 -independent. There is, however, no need to perform a 
modified pole model analysis to extract the 7r°-?7 matrix element, since we can obtain 
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it directly from chiral perturbation theory, at least for the region of q 2 for which the 
low energy expansion is valid.* 

Next we turn to the ^-dependence of m| 8 . If we consider the change in the matrix 
element in going from q 2 = m 2 to q 2 = -m 2 (very nearly equal to the range over which 
the p - uj matrix element is extrapolated to get from the e+e" experimental data to 
the scattering region) we find that there is a decrease of 17%. One of the advantages 
of a model calculation would be that the model would provide also the higher order 
q 2 dependence; here we have access only to the leading (linear-in-g 2 ) dependence, 
though this dependence has (apart from the uncertainty in the overall scale, a feature 
subject to future improvement) the advantage of theoretical reliability, in the sense 
that the structure of the low-energy expansion of chiral perturbation theory is directly 
related to QCD and its convergence to one-loop has been tested for a large number 
of processes 19 . We see that the q 2 dependence of the mixing in the meson propagator 
is certainly not negligible. While one would like to quote numbers on how much this 
affects, eg., the uncertainty in the n° - r\ contribution to the ^4 = 3 binding energy 
difference (in the conservative sense described above), this is not possible because 
the binding energy difference is, in fact, sensitive to spacelike q 2 values considerably 
outside the range of validity of the low-energy expansion 3 . In order to make further 

* It should be noted that the analysis of Refs 6a,6b which uses the quark model, 
together with certain results valid in the chiral limit, to evaluate, eg., the tt° - rj 
matrix element, contains non- vanishing EM contributions which are incorrect. The 
EM contribution to tt° - -q mixing actually vanishes in the chiral limit 13 . The error 
results from neglecting a class of graphs ("EM penguins of the second kind" 19 ) which 
are required to properly implement all of the EM chiral constraints of Ref 13, in the 
chiral limit. The numerical consequence of this error are, however, not particularly 
significant. 
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progress one should, therefore, use the chiral results embodied in Eqs (7) and (8) 
to constrain model calculations. If an extension of the calculation of Ref 11 to the 
pseudoscalar sector were to reproduce the ^-dependence of the ir°-ri matrix element, 
this would considerably enhance our confidence in what it has to say about the running 
of the p - u> matrix element. 
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Table 1. The constants a Q , 01 as a function of (Am K ) EM * 



(Am K ) EM (MeV) 


a (GeV 2 ) 


ai (GeV 2 ) 


1.3 


.00335 


.00028 


1.6 


.00354 


.00030 


1.9 


.00372 


.00032 


2.2 


.00390 


.00033 


2.5 


.00408 


.00035 


2.8 


.00426 


.00036 



* All notation as in the text. a D and a\ are as defined in Eq (8). 
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